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Abstract 

We compute energies and energy densities of static electromagnetic flux tubes in three and 
four spacetime dimensions. Our calculation uses scattering data from the potential induced by 
the flux tube and imposes standard perturbative renormalization conditions. The calculation is 
exact to one-loop order, with no additional approximation adopted. We embed the flux tube in 
a configuration with zero total flux so that we can fully apply standard results from scattering 
theory. We find that upon choosing the same on-shell renormalization conditions, the functional 
dependence of the energy and energy density on the parameters of the flux tube is very similar for 
three and four spacetime dimensions. We compare our exact results to those obtained from the 
derivative and perturbation expansion approximations, and find good agreement for appropriate 
parameters of the flux tube. This remedies some puzzles in the prior literature. 
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1 Introduction 



Flux tubes in QED coupled to fermions exhibit a number of interesting phenomena, such 
as the Aharonov-Bohm effect [T|, its consequences for fermion scattering [2 , parity anoma- 
lies fl], formation of a condensate j3j, and exotic quantum numbers El U\- The same 
(non-perturbative) features of the theory that give rise to these unusual phenomena make 
it more difficult to analyze this system with standard techniques, especially in calculations 
that require renormalization. The analysis in ref. jH] and the world line formalism in ref. 
have addressed some of these issues. Here we provide a comprehensive approach drawing on 
techniques from scattering theory to analyze the one-loop energy and charge of this system. 

Our primary motivation for this analysis is to shed light on vortices in more complicated 
field theories, especially the Z-string in the standard electroweak theory PU|- The Z-string 
is a vortex configuration carrying magnetic flux in the field of the Z-gauge boson. If a 
network of such objects existed at the electroweak phase transition, then it would be one key 
ingredient in a viable mechanism for electroweak baryogenesis without a first-order phase 
transition. Since the classical Z-string is known to be unstable jTTj, this scenario would 
require stabilization via quantum effects ^2], perhaps by trapping heavy quarks along the 
string. We also expect that by extending our results to the Abelian Higgs model, they 
could be applied to Abrikosov flux tubes in Type II superconductors [T3] or supersymmetric 
Abelian Higgs models [Tl] . 

We compare the one-loop energies and energy densities of electromagnetic vortices in 
D = 2 + 1 and D = 3 + 1 spacetime dimensions. The classical calculation is of course the 
same in the two cases. The quantum corrections to the energy could possibly be very different 
jH] because of the different divergence structure. In D = 3 + 1, the bare one-loop energy is 
divergent and only after we impose renormalization conditions do we get a finite result. In 
D = 2 + 1, in contrast, the bare energy is finite. However, a comparison between the two 
dimensionalities is sensible only when we use the same renormalization conditions, which 
requires a finite counterterm in the D = 2 + 1 case to keep the photon field normalization 
the same. Without this finite renormalization, the D = 2 + 1 and D — 3 + 1 energies are 
qualitatively different. But after proper renormalization, we find that both the energies and 
energy densities are closely related. 

We also study this problem to get a handle on several technical issues associated with the 
computation of the one-loop energy of a flux tube. As described in Section 13.21 an efficient 
way to compute the energy is to use scattering data of fermions in the background of the 
flux tube. However, vortex configurations give long-range potentials, which do not satisfy 
standard conditions in scattering theory [T3] , which would guarantee the analytic properties 
of scattering data. This leads to subtleties that are discussed in Section HJ We show that 
these puzzles arise only because an isolated flux tube is unphysical, and once a region of 
return flux is included, the scattering problem is well-defined and the puzzles disappear. In 
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the limit where the return flux is infinitely spread out, the energy density becomes entirely 
localized at the original flux tube. 

The paper is organized as follows: In the following Section we describe the theory and 
introduce electromagnetic flux tubes. We outline the computation of their classical and 
renormalized one-loop vacuum polarization energies in Section 3. We describe puzzles in this 
calculation that arise from an isolated flux tube in Section 4 and explain how these puzzles 
can be resolved by appropriate embedding of the flux tube in Section 5. We present numerical 
results for the energies, energy densities and charge densities in Section 6. In Section 7 we 
summarize our results and provide an outlook on related studies. Three Appendices give the 
technical details needed for the computation of quantum contributions to the energy and 
charge densities. 



2 The Theory 

We consider QED in spacetime dimensions D = 2 + 1 and D = 3 + 1, with a four-component 
fermion field, ip, in both cases. The Lagrangian density is 

C [D) = —F^F^ + i>(ifi+ e4-m)^j + , (1) 

where the indices fj, and v run from to D — 1 and CqQ is the counterterm Lagrangian in 
D spacetime dimensions. In D — 3 + 1, this Lagrangian describes electromagnetism with 
a single four- component fermion of mass m and charge e. In D — 2 + 1, we have parity- 
invariant electromagnetism with two flavors of two-component fermions of equal mass m and 
charge e. We are interested in static magnetic flux tubes. These are localized, cylindrically 
symmetric magnetic fields (pointing in the z direction in D = 3 + 1), with a net flux F 
through the xy-plane. In the radial gauge, flux tubes arise from vortex configurations of 
gauge fields: 

A) = 0, A=£-f(r)e ip , (2) 

where r 2 = x 2 + y 2 measures the planar distance from the center of the vortex. The radial 
function /(r) goes from Otol. In .D = 2 + 1 the magnetic field is the radial function 

w 2vrr dr ' w 

while in D = 3 + 1 it is the vector field B (r) = B(r)e z . For small r, we must have f(r) 
going to zero at least quadratically for B(r) and A(r) to be non-singular. 

We often find it convenient to specify the flux in units of 27r/e and define the dimensionless 
quantity 
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In most cases we will take the example of a Gaufiian flux tube of width w, 

f G (r) = 1 - e"^ 2 so that B G (r) = 5 G (0) e -^ 2 , (5) 

whose flux is 



F G = ™ 2 5 G (0). (6) 



3 The Energy 



We would like to to compute the total energy £(3) for D = 2 + 1 and the energy per unit 
length for D = 3 + 1, to one-loop order. In general, we would have to evaluate the 
appropriate matrix element of the energy momentum tensor. However, as we review in 
Appendix [Bl for static configurations the total energy is simply given by the negative one- 
loop effective action per unit time. Since the theory is Abelian, the photons do not have 
self-interactions and the one-loop effective action is obtained by integrating out the fermion 
field. The photon fluctuations start contributing only at two loops and higher, which we 
ignore. 

To one loop order, the total energy is the sum of the classical energy, E c \, and the fermion 
vacuum polarization energy, E^ , 

= E cX + E® . (7) 
For both cases, D = 2 + 1 and D — 3 + 1, the classical contribution is 



E cl = i J SrB\r) . 



We extract the renormalized fermion vacuum polarization energy from the one-loop effective 
action 

E®. = lim - [ In Det(i^+ efi-m)- In Det(^- m)] + E§^ , 

E^= lim [lnDet(^+e4-m) - lnDet(^-m)] + , (9) 

by dividing out the arbitrarily long time interval in both cases, and the arbitrary length of 
the vortex in D = 3 + 1 dimensions. Note that the energy is always defined relative to a 
background where the electromagnetic fields are zero everywhere. 
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3.1 Renormalization 



The counterterms introduced in eq. originate from the standard counterterm Lagrangian, 
£ct — ~^i£-Fuv Since we do not consider electric fields, the counterterms are quadratic in 
the magnetic field 

Effl — Jfx*. (10) 

To fix the counterterm coefficient, C^ D \ we choose on-shell renormalization conditions, which 
require that the residue at the q 2 = pole of the photon propagator is unity to ensure that 
one-photon states remain normalized. 1 We obtain 

c m r(2-f) 

3(4vr) D / 2 m*- D 1 ' 



which becomes 



C {3) = ~^—, (12) 
bum 

in D = 2 + 1 dimensions, and the dimensionally regularized quantity 

oM --5.(!-' + '»S)- (13) 

for dimensions approaching D = 3 + 1 , where e = 4 — D . 

The counterterm, £ ct renormalizes the bare photon field and fermion charge, 

A» = (l + CM)- 1 t 2 A£ m and e=(l + C^) 1 / 2 e bare . (14) 

In the absence of photon fluctuations, the bare fermion field and mass do not get renormal- 
ized. In D = 3 + 1 the counterterm is divergent and combines with the divergent functional 
determinant to give a finite, renormalized one- loop energy. In D = 2 + 1, the counterterm 
is finite, but it still must be taken into account to ensure that the photon field remains 
canonically normalized. Since we want to compare the energies in D = 2 + 1 and D = 3 + 1, 
we must impose identical renormalization conditions in the two cases. 



3.2 Functional Determinant 

We use the phase shift approach [TBI EI to compute the functional determinant in eq. (JSJ) 
exactly. The fermion vacuum energy is given by the renormalized sum over the shift in the 

1 Although fcrmions are logarithmically confined in D = 2 + 1, asymptotic photon states can still be 
observed experimentally through photon-photon scattering. 
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zero-point energies of the fermion modes due to the presence of the background magnetic 
field. This calculation comprises a sum over bound state energies and an integral over 
the continuum energies weighted by the change in the density of states. These quantities 
are obtained from the Dirac equation for the fermion fields in the background of the flux 
tube. The explicit form of this equation is given in Appendix |X1 eqs. (JoTij) and (pTTj). Since 
the configuration that characterizes the flux tube is cylindrically symmetric, the fermion 
scattering matrix may be decomposed into partial waves labeled by the z-component of the 
total angular momentum, M. Then, in D — 2 + 1, the change in the density of continuum 
states is given in terms of the phase shifts of the fermion scattering wave functions in the 
flux tube background, 

A PW = E - 1 [<WW] , (15) 



7r dk 



where k is the magnitude of the momentum and a refers to all other discrete labels. For a 
given momentum k there are two energy eigenvalues ±^/k 2 + m? and two spin states. This 
four-channel scattering problem can easily be diagonalized, yielding four identical phase shifts 
(with the exception of the threshold at k = 0, as discussed below). Therefore summing over 
the discrete labels merely gives a factor of four. We denote the result of this sum as Suik). 
In D = 3 + 1, we must also integrate over momenta in the z-direction [IB] . The details of 
the phase shift calculation are given in Appendix 1X1 

To renormalize the continuum integral, we subtract the first two terms in the Born expan- 
sion of the phase shift and add back in the corresponding energy in the form of the two-point 
Feynman diagram (a fermion loop with two insertions of the background potential). The 
details of this calculation are described in Appendices IA1 and IH1 In D — 3 + 1, the Born 
subtraction renders the integral over the continuum energies convergent. The ultraviolet di- 
vergences are isolated in the Feynman diagram, which when combined with the counterterm, 
eq. ([TU]) . gives the properly renormalized finite result. In D = 2 + 1 the Born subtraction is 
finite, but implements the on-shell renormalization condition. 

Our final expression for the renormalized fermion vacuum energy is the sum of the sub- 
tracted contribution of the bound and continuum modes, the Feynman diagram, and the 
counterterm: 

E^=E^+EiS+E^. (16) 

The contribution from the Born-subtracted mode sum can be expressed through scattering 
data, 



*f - -h E {?, '4 + ^ - 4) - h [ dkk ^ £ • < 18 > 

j \ / JU M 
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where Uj are the energy eigenvalues of the discrete bound states. (Even though the con- 
figurations that we consider here do not have such bound states, we have included their 
contribution for completeness.) We have defined 



5 M {k) = 5 M {k)-5 M \k)-8^{k) 



f( 2 ) ( 



(19) 



where 8$ denotes the i th term in the Born series expansion of the exact phase shift, 8m- 
Since we have subtracted the two leading contributions of the Born series, the remainder, 
E$ D ^ comprises the sum of the third and higher order pieces. 

Next we add back the two leading Born terms in form of Feynman diagrams. The first or- 
der diagram vanishes by Furry's theorem. In dimensional regularization the unrenormalized 
second-order Feynman diagram energy is 

df(r) 



E 



(D) 
FD 



(4tt) d / 2 



dp 



dr- 



dr 



-J (pr) 



dx 



x(l - x)pT(2 - D/2) 
[m 2 + p 2 x(l — x)] 2 ~ D l 2 



(20) 



Combining this expression with the counterterm energy, eq. (jlUJl . gives the renormalized 
Feynman diagram energy at the physical spacetime dimension. The result is to impose 
perturbative renormalization conditions on this non-perturbative calculation. 



4 Subtleties of Configurations With Net Flux 

We have described the general computation of the renormalized one-loop energies of mag- 
netic flux tubes. Before proceeding with the calculation, however, we must address subtleties 
in the calculation of the phase shifts at zero momentum. These questions arise because the 
background potential does not satisfy standard asymptotic conditions in scattering theory. 
As is well-known from the Aharonov-Bohm effect, even if the magnetic fields are exponen- 
tially localized, the vector potential falls only as 1/r, which is too slowly for many of the 
standard theorems of scattering theory to hold [T3]. As a result, the phase shifts are discon- 
tinuous at threshold. Although this discontinuity does not cause problems for calculating 
the vacuum polarization energy, because the integrands in eqs. (|17|) and (J 18)) remain smooth, 
we can no longer determine the number of bound states through Levinson's theorem ^H] ■ 

As described in Appendix |X1 we calculate the phase shifts in the background of a flux 
tube by deriving second-order differential equations from the Dirac equation. The asymptotic 
behavior as r — > oo of the coefficient functions in these differential equations is extracted 
from eqs. (|53|) in with h(r) = JF and i — M — |, 



. m . (r) + g^0 + ^_(^LtizfV) 9 « (r) = , (2i) 
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where primes denote derivatives with respect to r. For large r we do not have free Bessel 
differential equations whose index equals the angular momentum; rather the equations de- 
scribe an ideal flux tube with profile function f(r) = 1. If extended to the origin, such an 
ideal flux tube would generate a singular background potential. 

Since the centrifugal barrier is shifted by the amount of the flux, the regular solutions to 
the asymptotic differential equations (|21j) are Bessel functions with a correspondingly shifted 
index, 

9 { Pij) = J\ M -i-F\(kr) and gf\r) = J {M+ i_ rl (kr) . (22) 

From the asymptotic behavior of the Bessel functions, we read off a phase shift relative to 
the trivial configuration = [201 

AS M {k) = 7T (\M - i| + \M + i| - \M - i - T\ - \M + ^ - F\) , (23) 

where we have summed the phase shifts associated with the two radial functions, which 
corresponds to summing over the spin degrees of freedom. For a radially varying profile 
function /(r) with lim r ^ 00 /(r) = 1, the second-order equations (JH3|) allow us to compute 
phase shifts relative to the ideal flux tube, 2 <5Af,ideai, to which ASm must be added to obtain 
the phase shift relative to the trivial configuration: 

8m = <Wdeal + A<5 M . (24) 

However, for k = 0, there is no r large enough where the asymptotic form of the Bessel 
functions holds. Hence equation (|23|) is valid only for k > and and it is unclear what 
A5jvr(0), and hence £m(0), should be. 3 

We also note that the second-order differential equations are not equivalent to the Dirac 
equations for k = 0. For uo — ■> +m the lower component of the spinor becomes identically 
zero, while for uj — > —m the upper component is zero, 4 so the elimination of one component 
in favor of the other becomes singular. For most applications, this is not a problem because 
the phase shifts have a smooth limit as k — > 0. However, as discussed in Appendix IA.21 
the long-range potentials associated with flux tube configurations allow for bound states 
with zero binding energy, which introduce discontinuities in the phase shift at threshold. 
This singularity is reflected in our numerical analysis of the phase shifts. For k significantly 
greater than 0, the phase shifts that we obtain from the two radial functions gf\r) and 

2 The corresponding scattering wave-functions approach H^ ±3 __ :F {kr) — e 2lSM - idoal H^ 1 > ±L _ :F (kr) at large 

radii, where are Hankel functions. 

3 See refs. |22 for discussions of the relation between this discontinuity of the phase shifts and the 
anomaly in D = 2 + 1 . 

4 We use the Bjorken-Drell convention for Dirac matrices. 
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gf\r) are identical, as they should be. In the limit k — * 0, we find that for M > | only the 
differential equation for gf\r) converges to the phase shift, while the equation for gp gives 
discontinuities in Suik) because of the threshold bound states. The discontinuity is smeared 
out from k = because in practice we integrate the second-order equations starting from a 
large-distance cutoff, rather than infinity. For M < — |, situation is the same, with the roles 

of gp and gf' reversed. For small but nonzero k, the phase shifts (without singularities) 
converge to 

lim^ f m-vrl |M + ^" |M+ ^-- F| fOT (25) 
hmd M (k)-7T^ M _ i|_| M _|_ JC -| for M <-\- {2b) 

According to ref. ^Hj we should be able to compute the number of bound states from 
the the phase shifts at k = 0, 

Nm = ~ [M w +»n) + M"> -> -m)] = - lim M*0 (26) 

7T 7T fc^O 

since lim^oo S^ik) = 0. With the help of an explicit example we will now show that this 
result does not hold in the present case. For definiteness we consider T = —1.2. From 
eq. (jHHj) in Appendix IA.2I we conclude that there is a threshold state at u> = —m with 
M = -1/2. However eq. flU yields 

1.6 for M = -\ 
Nfj = { 2.4 for M < -\ (27) 
-2.4 for M > i 



5) 



and we have a discrepancy between the number of states that leave the continuum (jV 
and the number of bound states (1). 

These problems can be seen more concretely by considering the case of D = 2 + 1 with 
a single two-component fermion. Then the net fermion charge in the background of the flux 
tube is —J r /2, which is connected to the parity anomaly jjj. Thus charge conservation is 
violated as the flux tube is turned on. Furthermore, if we hold the flux fixed but spread 
the magnetic field out over a larger and larger region the total energy of the flux tube 
approaches zero (as shown in Section ^. 1.2j) . leading to the implausible conclusion that there 
exist charged, massless objects in the theory and the fundamental fermion cannot be stable. 

The puzzles addressed in this Section do not necessarily prevent us from calculating 
the vacuum polarization energy, because the small momentum region is suppressed in the 
integrals in eqs. (|171 fTHj) and we have methods other than eq. (|26[) available to find the 
number of bound states and compute their binding energies. However, to be certain of the 
consistency of our results, we turn to a method to eliminate these puzzles by extending the 
background potential to satisfy the standard conditions of scattering theory. This approach 
will then also allow us to directly generalize the methods of ref. to compute energy 
densities. 
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5 Embedding 



All the issues discussed in the previous Section can be resolved by recalling that it is not 
possible to create a configuration carrying non-zero net flux starting from nothing. In 3+1 
dimensions we have the Bianchi identity 

6^8^ = 0. (28) 

It is a mathematical identity and not an equation of motion derived from the action principle. 
In terms of electric and magnetic fields, it gives us two of the Maxwell equations: 

V-B = 0, — = -VxE. (29) 
dt v ; 

The fact that the magnetic field is divergenceless requires all magnetic flux lines to be closed 
and there can be no net flux through a plane unless closure of the flux lines is arranged at 
spatial infinity, which may not be appropriate to compute the vacuum polarization energy 
since that requires integrating over all space. In 2+1 dimensions, only the second of the two 
equations remains in the form 

OB 

— = -d x E y + d y E x . (30) 

Integrating this equation over space, it is clear that the magnetic flux is time independent 
for localized fields and it is not possible to create a net flux starting from no flux. 

Therefore, we consider flux tube configurations with a return flux such that all flux lines 
are closed and the total flux is zero. We say that the flux tube is embedded in a physical 
no net flux configuration. Any apparently missing states or charge can be accounted for 
by considering the corresponding quantities localized around the return flux. Once the flux 
tube is embedded in a no net flux configuration, the resulting potentials in the second order 
differential equations for gp and gp fall into the class of potentials that are well understood 
in scattering theory. (Ref. [S] uses similar ideas as part of a complementary approach to the 
scattering problem.) The separation between the locations of the flux tube and the return 
flux should be large compared to the extent of both the flux tube and the return flux. As 
we will demonstrate, the vacuum polarization energy of the configuration with zero net flux 
approaches a well-defined limit as the return flux is sufficiently separated from the flux tube 
and is sufficiently diffuse. This limit represents the energy of the flux tube alone, which we 
can verify by computing the corresponding density over the volume of the flux tube. We find 
that this result agrees with the energy obtained by using the phase shifts in the background 
of an isolated flux tube, ignoring the subtleties discussed in the previous Section. 

We now describe the above embedding in more detail. We consider a special subset of 
return flux configurations with magnetic field Br for which the extension, wr, is proportional 
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to the distance, R, from the flux tube. Since then Br is characterized by a single length scale, 
and the flux is held fixed, simple scaling arguments based on the perturbative expansion show 
(see Section loM.2|) that the classical energy goes to zero like 1/R 2 and the one-loop energy 
with on-shell renormalization goes like 1/R A in both D = 2 + 1 and D — 3 + 1. Thus the 
energy of the no net flux configuration should approach the energy of the flux tube alone 
as R — > oo. To show that this value corresponds to computing the energy using the phase 
shifts in the problem with no return flux, we consider the Gaufiian flux tube, Bq(t), defined 
in eq. (J5J) with flux Fa- We take 

B R (r) = -, , (31) 

irR 2 (1 + 256 (r 2 /R 2 - if) (vr/2 + arctan(16)) 

which originates from the profile function 

. . arctan [16(r 2 /R 2 — 1)] + arctan(16) , . 

Mr) = vr/2 + arctan(16) (32) 

for our return flux. By construction /ij(r) — > — 1 as r — > oo such that /G( r ) + /n( r ) — > in the 
same limit, with /b( r ) defined in eq. (|HJ). Since this sum also vanishes as r — > 0, it generates 
a scattering potential satisfying standard conditions in scattering theory. Analogously for 
the magnitic field, we have the no net flux embedding 

B (r)=B G (r)+B R (r). (33) 

The total energy consists of three parts: the classical energy, the renormalized Feynman 
diagram energy and the phase shift energy. We consider these different energy contributions 
for B as a function of R. We choose Tq = 4.8 and w = 1/m to fix the flux tube, Bq. Using 
the expression for the classical energy, eq. (JBJ), it is straightforward to verify that 

lim E cl [B ]=E d [B G ]. (34) 

R^oo 

In Fig.^ we show the renormalized Feynman diagram contribution to the energy , E^+Eqt 
as a function of R as computed from eqs. ()20j) and (jlOj) for the no net flux configuration, 
eq. (J3l?|) . In both cases, D = 2 + 1 and D — 3 + 1, we find that the limit is well saturated 
for R > 10/m. In Fig. |21 we show the integrand of the phase shift part of the energy, 
eqs. (fTTl ITSj) . in both the embedded configuration and the flux tube configuration, for two 
fixed values of R. The details of the calculation of the phase shifts of the non-zero flux 
configuration may be found in Appendix IA.lt We observe that the integrands disagree only 
for small values of the momentum, at which the fermion is sensitive to the presence of the 
spread-out return flux. The integrand for the embedded configuration oscillates around the 
integrand of the flux tube configuration with an amplitude that decreases as k increases. The 
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Figure 1: Renormalized D = 2 + 1 and D = 3 + 1 Feynman Diagram energies in appropriate 
units of m as functions of the return flux radius for w = 1/m and T = 4.8. The solid lines 
correspond to the energies without the return flux. 

region of disagreement gets pushed to smaller values of k as R increases. In the limit R — > oo, 
we find that the phase shift contribution to the energy is identical for the embedded and the 
flux tube problems. Analogous results hold in D = 3 + 1 because the integrand differs from 
that in D = 2 + 1 only by a background- independent function of k, cf. eqs. (fT7|) and (fTSj) . 

Since the energy in the embedded problem approaches the energy in the isolated flux 
tube problem, for computational efficiency we can use the phase shifts in the background of 
the isolated flux tube to compute the energy, with the implicit understanding that there is 
a spread-out return flux present that does not affect the energy. 



In this Section we discuss our numerical results in more detail. The central issue will be the 
energy of the flux tube. We will also consider its energy and charge densities, as obtained 
from the embedded problem, and compare our results to the derivative expansion. 



6 Results 
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Figure 2: Integrand of the phase shift contribution to the energy in D = 2 + 1 for R = Q/m 
(left panel) and R = 26/m (right panel). The solid lines correspond to the isolated flux tube 
problem and the dashed lines to the embedded problem. 

6.1 Energy 

As already argued in the previous Section, we can compute the energy of a magnetic flux 
tube in three and four dimensions either independently or by embedding it in a zero-net-flux 
configuration. We then use scattering data to compute the vacuum polarization energy. Of 
course, there are other (approximate) methods to compute the vacuum polarization energy 
of such a configuration, most notably the derivative expansion [1J|^J|22], which represents 
an expansion in derivatives of the magnetic field. However, this expansion has recently been 
criticized, in particular in the four- dimensional case jUJ- Since our approach is exact, we can 
easily test the validity of the derivative expansion. In addition, we may compare our results 
to a perturbative expansion in the magnitude of the magnetic field. For the Gaufiian flux 
tube, eq. (jHJ) we can formulate both expansions by varying the width, w under appropriate 
constraints. 

6.1.1 Derivative Expansion 

As its width approaches infinity, a Gaufiian flux tube with a fixed value of magnetic field at 
its center approaches a spatially constant magnetic field, for which the derivative expansion 
is reliable. Therefore, to compare with the derivative expansion, we compute the total energy 
of configurations with varying w and fixed -Bg(O) for the Gaufiian magnetic field, eq. (JSJ). 
The flux, eq. (jUJ), is not held fixed in this approach. 

In D = 2 + 1, the derivative expansion of the unrenormalized one- loop vacuum polariza- 
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tion energy is given by |U 121] 



77.(3) _ 77.(3) _ 77.(3) 77.(3) fOC^ 

to next-to-leading order in the derivative, where 

4to = / **^T [ dse-™ 2 ^s-^ (coth( S )-±) , 

4e,2 = \J ^[V(eE)l 2 |47reEr 3 / 2 |^d S e-- 2 /' eg l g -^ ^ggM , (36) 

Note that the above expressions are for four-component fermions in the loop, which give twice 
the result for two- component spinors. Scaling the spatial integration variable by £ = x/w 

straightforwardly yields -E^eo ^ u ' 2 anc ^ ^de2 0(1 ' u;0 - The counterterm contribution to the 
energy is also proportional to w 2 . We stress that keeping B(0) fixed is important to obtain 
this simple power law behavior. (Taking the limit differently, for example with the flux fixed, 
the proper-time integrals in eqs. (j3Tj) would induce more complicated dependences onto.) 
For the special case of the Gaufiian flux tube we find 

{3) _ e 2 B 2 G (0)w 2 m] 
We add the counterterm to get the renormalized derivative expansion 

E llc,DE = E CT + E< DE,0 + E DE,2 + ■■■ ■ (38) 

It is straightforward to verify that this expression satisfies the on-shell renormalization con- 
dition that we imposed on our scattering data result in subsection 13.11 The first two terms 
on the right hand side of eq. (|3*K|) are proportional to w 2 , while the last term is independent 
of w. All omitted terms contain more derivatives, and thus they vanish as w — > oo. 

In D = 3 + 1, the next-to-leading order derivative expansion of the renormalized one-loop 
energy is 

J vac,DE "^DE^ 



77,(4) 77.(4) 77.(4) ( „Q\ 



where j2Hj 

4to = / 1°° dse ' sm 2/leBls ' 2 ( coth (*) - \ - 1) , 

E nl,2 = - / rf 2 x|V(efi)| 2 |327T 2 e5r 1 (40) 
x J™ dse- sm2 '\ eB \ ^l-4coth 2 S + 3coth 4 S + ^^(l-coth 2 S )^ . 
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Figure 3: Renormalized one-loop energies in D = 2 + 1 (left panel) and D = 3 + 1 (right 
panel), for fixed values of the magnetic field at the origin, as a function of the width of the 
GauBian flux tube. The lines correspond to the derivative expansion to lowest nontrivial 
order in the derivative. From top to bottom, eSo(O) / m 2 = 1.1,2,2.5,3,4,5. 

Again we have oc w 2 and E^ 2 oc w°, with higher orders going as inverse powers of w. 

In Fig. El we compare our exact result E^} with the derivative expansion approximation 
^iac de f° r different values of ei?G(0), and find excellent agreement for widths larger than 
1. There appears to be no qualitative difference D = 2 + 1 and D — 3 + 1. However, had 
we not included the counterterm in D — 2 + 1, as is done in jU], we would be effectively 
imposing a different renormalization condition from the one in D = 3 + 1, and would see a 
qualitative difference between these two cases. As we will see below, the extra contribution 
from the off-shell scheme can dominate the underlying result. In particular, if this scheme is 
consistently imposed both exactly and in the derivative expansion, this dominant contribu- 
tion can obscure differences between the exact result and the derivative expansion in cases 
where the derivative expansion is not valid. 

6.1.2 Fixed Flux 

We next turn to such a limit, which is the case of configurations for which the flux, rather 
than the magnetic field at the origin, is held fixed as the width approaches infinity. This 
limit was considered in ref. [§]. In this case the energy of the flux tube goes to zero as w 
goes to infinity, and it is the perturbation expansion, rather than the derivative expansion, 
that becomes exact. 
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First consider the classical energy, eq. (JSJ), 



e ci =^ r~ ( dfG{r) y , (4i) 

e 2 J r \ dr J e 2 w 2 ' 

for the GauBian flux tube, eq. (0). For large widths, it goes to zero like 1/w 2 . In this limit, 
the magnetic field becomes weak, as can easily be seen from eq. ©• Hence the dominant 
contribution to the vacuum polarization energy comes from the two-point function. In 
.0 = 2 + 1, the unrenormalized two-point energy can be expressed as a series in 1/w 2 . The 
leading term proportional to 1/w 2 turns out to be exactly equal to minus the counterterm 
energy in the on-shell subtraction scheme. Thus only the subleading term survives in the 
renormalized two-point energy. For a GauBian flux tube, 

bmw z \ bm^w 1 
Eg + Eg = -^r^ + 0(l/w e ) . (42) 

For large widths, renormalization not only changes the sign of the one-loop energy, but also 
the rate at which zero is approached. The renormalized results in D = 3 + 1 are similar to 
the results in D = 2 + 1: 



£ ™ + £ o--30^ + ° (1/ra6) (43) 



for the GauBian flux tube. In the limit w — > oo, the total energy for both D = 2 + 1 and 
O = 3 + 1 vanishes. This result enables us to introduce a return flux in the embedded 
problem without changing the energy. 

While such fixed flux configurations are well approximated by the perturbative expansion, 
they are inappropriate to test the derivative expansion, as was done in ref. 0. Since the 
magnitude of the magnetic field at the center of the vortex then depends on the width, so 
too will the proper-time integrals in eqs. (j36ll4*Uj) . and there is no longer any reason to expect 
that the the two-derivative contribution is less important than the no-derivative contribution 
for large widths. 

In Fig. El we plot the exact one-loop energies for various values of the flux as a function 
of the width. We normalize the energies in units of T 2 so that all differences between the 
different fluxes are due to three-point and higher contributions. We compare the energies 
with the leading order two-point energies and find good agreement for large widths. As the 
flux increases, we need to go to larger widths to get a weak magnetic field everywhere, cf. 
eq. Hence the agreement with the leading order two-point function contribution to the 
energy sets in at larger widths when the flux increases. 
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Figure 4: Renormalized fermion vacuum polarization energy in units of J 1 " 2 as a function of 
the width, for various fixed values of the flux JF (2.5, 4.5, 6.5, 8.5 from top to bottom) in 
the GauBian flux tube. The full line represents the two-point function contribution. The left 
panel is for D = 2 + 1 and the right panel for D = 3 + 1 . 



6.1.3 Further Comments 

Methods similar to our phase shift calculation have been used in ref. j2H|. However, in this 
calculation infinite quantities are used without proper renormalization, such as eq. (26) of 
j2bj for D = 3 + 1. Although convergent results are reported, we believe that they are 
an artifact of taking a low, fixed upper limit in the sum over channels. Instead, this limit 
should scale with the momentum k. By keeping more terms in the sum, one should see this 
divergence. 

n Ref. jHj a step function background B(r) was studied in D = 2 + 1 dimensions. In 
order to separate divergent contributions to the vacuum polarization energy the authors 
of Ref. jH] extracted the asymptotic behavior of the Jost function and identified that with 
pieces in the heat kernel expansion. Instead of using on-shell renormalization conditions, the 
renormalization is defined by requiring that the renormalized vacuum polarization energy 
vanishes in the limit where the fermion becomes infinitely heavy. However, this requirement 
is not sufficient to uniquely determine the renormalization scheme, as has been noted in |2*7] . 
As a result, we cannot make a quantitative comparison with our calculation, but qualitatively, 
our results seem to agree in the sense that they also find a negative vacuum polarization 
energy. Since different renormalization conditions amout to different contributions from the 
positive definite counterterm, eq. (jlUj) . this agreement may be not be conclusive. 

6.2 Energy Density 

So far we have only discussed the total vacuum polarization energy of the flux tube. We have 
argued that it is preferable to consider a flux tube embedded in system with no net flux, 
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Figure 5: The vacuum polarization energy densities of configurations with zero net flux as 
functions of the separation R = 8, 10, 12/m between Bq and Br. We choose eB(0) = Am 2 
and w = 1.5 /m. The left and right panels show D = 2 + 1 and D — 3 + 1, respectively. 



where the region carrying the return flux is well separated from the flux tube core. In this 
Section we make this argument more precise by studying the radial density of the vacuum 
polarization energy. We will observe that the contribution to the energy density from the 
central flux tube and the return flux can clearly be distinguished, giving a unique definition 
for the energy of the central flux tube alone. Numerical evaluation shows that it equals the 
energy of the non-embedded flux tube that we computed in Section El We will also take 
the opportunity to compare our results for the energy density with the derivative expansion 
approximation, which has recently been discussed in ref. j^]. 

We define the radial energy density in terms of the renormalized vacuum expectation value 
of the energy-momentum tensor T^ix) by e(r) = 27rr(Too(x)). In Appendix iBl we discuss 
Tqo(x) in more detail, with particular attention to total derivative terms. The radial energy 
density e(r) has been normalized so that the total energy is = dre(r). In figure El 
we show the energy density for the configuration defined in eq. (J33|) . As the separation R 
between Bq and Br increases, the contribution from the return flux region decreases, going 
to zero as R — > oo. In addition, we observe that the energy density localized around the flux 
tube at r = remains unchanged when we increase the separation. That part of the energy 
density does not depend on whether the return flux configuration is included once Bq and 
Br are well separated. Finally, the intermediate region in which the density e(r) vanishes 
becomes more pronounced. Once the region with e(r) w is large enough we can clearly 
distinguish between the energy densities due to B G and B R . 

In figure |H1 we show the renormalized Feynman diagram and scattering data contributions 
to the energy density. We observe that for small widths of the central Gaufiian flux tube, the 
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Figure 6: Feynman diagram (dashed) and scattering data (dotted) contributions to the 
energy density for conhgurations with zero net flux. Their sum (full) is the full vacuum 
polarization energy density. Two cases are considered, w = 0.75/m (left panel) and w = 
1.5/m (right panel) , both with eB(0) = Am 2 and R = 8/m, for D = 2 + 1 (top) and D = 3 + 1 
(bottom). 
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Figure 7: Renormalized Feynman diagram contributions to the energy density for configu- 
rations with and without return flux for w = 1.5 /m, eB(0) = 4m 2 and R = 8/m, in the case 
of zero net flux. The left and right panels show D = 2 + 1 and D = 3 + 1 respectively. 



scattering data contributions are negligible, while for larger widths they contribute as much 
as 30% to the total energy. Note that B(0) is kept fixed. We also see that the scattering 
data contribution to the energy density essentially vanishes in the return flux region. This 
numerical result reflects a cancellation when summing over many orbital angular momentum 
channels; the individual channels themselves contain sizable contributions to the energy 
density in the vicinity of r = R. This result is to be expected: Since the return flux region 
extends to large radii, large angular momentum channels contribute to the energy density 
there. We find that several hundred channels must be summed for convergence. 

The computation of the Feynman diagram contribution to the energy density does not 
rely on any conditions for the background field. Since the energy densities are well approx- 
imated by their Feynman diagram contributions for small widths, we can compare energy 
densities for configurations with and without Br in cases with small w. This comparison is 
shown in figured 

Integrating the energy density over the region of return flux, we find that even at moderate 
separation it only contributes a small amount to the vacuum polarization energy. However, 
as figure illustrates, this result arises from the cancellation of significant positive and 
negative densities in this region. It is not surprising that the energy density of the return 
flux is nontrivial: We keep B(0) fixed, so as the width of the flux tube increases, its flux 
does too, according to eq. (jUJ). To have zero net flux, the amplitude of the return flux piece 
must increase when w gets larger and R remains unchanged, as can be seen from eq. (J31j) . 

To compare the results for the exact vacuum polarization energy density with the deriva- 
tive expansion, we must identify the energy density in that approximation. The derivative 
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expansion for the vacuum polarization energy, eqs. (j36l4U|) . originates from the expansion 
for the action. Hence simply omitting the radial integrals in those expressions can only be 
expected to yield an approximation to the action density. Formally we may write 

E^ = fdrr(^Je\B(r)), (44) 

with fi D \B(r)) to be read off eqs. (|3b14U|) for n = 0,2. As discussed in Appendix IB1 the 
action and energy densities differ by total derivatives, which vanish when integrated over 
space. These total derivates only affect the orders n > 2 because for constant magnetic 
fields both the energy and action density should be constant. Thus only E^ 2 is affected 
at next-to-leading order. The direct determination of this total derivate term involves all 
orders in perturbation theory and seems difficult to accomplish. We therefore introduce the 
parameter £ to define the energy density at next-to-leading order in the derivative expansion, 



43(r) 



/o (D) (^)) + 



dB(r) 
dr 



f { 2 D \B{r)) 



dr 



dB ( r ) AD) (U( \ \ 
r— — f y 2 >{B{r)) 



dr 



(45) 



where the case £ = corresponds to the the derivative expansion for the action density. 
Then we can fit for £ by comparing to the exact result e(r). 

For the derivative expansion to be applicable we need to consider increasing widths w 
with B(0) unchanged. As discussed above, for such configurations we may not omit the 
scattering data contribution and therefore we consider configurations with zero net flux. 
The comparison of the energy density is displayed in figure |H1 We find reasonable agreement 
with the derivative expansion for sizable w in the region of the central flux tube for both 
D = 2 + 1 and D — 3 + 1, where the magnetic field is slowly varying. Since we keep the 
position and the width of the return flux unchanged as we increase the width of the flux 
tube, we do not expect the derivative expansion to match the exact result in the vicinity of 
the return flux. From figure |H1 we observe that £ = —2 appears optimal for both D = 2 + 1 
and D — 3 + 1. This result confirms the assertion of Appendix El that the energy and action 
densities differ by total derivatives. 

In fig. El we display the normalized difference [e(r) — eDE( r )] /N with the normalization 
factor N = max{|e(r)|, r < 2w} for various widths w. Of course, we have to ensure that there 
is no overlap with the return flux when we consider large values of w. For w = 2.5/m (the 
largest value we consider) this requirement is well satisfied for R = 10/m. As expected, this 
normalized difference decreases as w increases. A non-zero difference extends to larger r as 
w increases, reflecting the growing extension of the central flux tube. Again, we observe that 
the inclusion of the total derivative term improves the agreement with derivative expansion 
approximation. As can be clearly seen from fig. El even though the difference between the 
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Figure 8: Comparison between the computed energy density for configurations with zero 
net flux and the derivative expansion approximation. We take w = 2.5/m, eB(0) = 4m 2 
and R = 10/m and consider the case without (£ = 0) and with (£ = —2) total derivative 
term, (|45|) . The left and right panels show D = 2 + 1 and D = 3 + 1, respectively. 



derivative expansion and the exact result for the density is quite sizable at small widths, it is 
close to a total derivative, which explains the good convergence of the derivative expansion 
for the total energy, even at moderate widths w. 

To summarize, we find reasonable agreement for the energy density between our exact 
results and the derivative expansion approximation at next-to-leading order for both D = 
2 + 1 and D = 3 + 1. Our result for the D = 3 + 1 case contradicts the findings of ref. [D]. That 
work uses the world-line formalism as a non-perturbative computation of the action density 
and finds a discrepency with the derivative expansion approximation in D = 3 + 1. We 
note that this discrepency cannot originate from the total derivative terms in the derivative 
expansion, since that work also considered the action density. Also, we have seen that the 
total derivative terms affect D = 2 + 1 and D = 3 + 1 similarly. 

On the other hand, in ref. a fixed flux sequence of magnetic fields has been employed 
to test the derivative expansion. We have already argued that a fixed peak magnetic field 
configuration should be used instead. Thus one might wonder why ref. [§] finds a good 
agreement for a fixed flux configuration in the D = 2 + 1 case. We ascribe this result to 
the choice of renormalization scheme, in which the finite counterterm has been omitted. For 
both the energy density and the total energy, we find that the inclusion of the counterterm 
leads to considerable changes in the case of D = 2 + 1. Since the counterterm only affects the 
second-order contribution to the energy and energy density, it is sufficient to consider the 
corresponding Feynman diagrams to study the influence of renormalization for D = 2 + 1. 
This comparison is shown in figure The renormalized energy density is about an order 
of magnitude smaller and of opposite sign than the unrenormalized one. Thus if one does 
not add a counterterm in D = 2 + 1, the leading order term £de,o contains a large, though 
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Figure 9: The difference between the exact result and derivative expansion for the energy 
density in the vicinity of the central flux tube as a function of its width. These differences 
are always normalized to maximal absolute value of the exact energy density in the flux tube 
region. The left and right panels show D = 2 + 1 and D = 3 + 1, respectively. Also we 
consider £ = (top) and £ = — 2 (bottom). Note the different scale of the latter. 
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Figure 10: Comparison of renormalized and unrenormalized 2 order Feynman diagram 
contribution to the energy density for D — 2 + 1. For simplicity we set Br = and take 
eB(0) = Am 2 and w = 1.5/m. 

unphysical, contribution to the energy and energy density. But then the comparision with 
the derivative expansion merely tests the agreement in this term, which masks the under- 
lying differences between the exact and approximate results. We thus conjecture that the 
authors of ref. jU] would also have found disagreement between their results and the deriva- 
tive expansion approximation in D = 2 + 1 if they had imposed the on-shell renormalization 
condition for their fixed flux background configuration. In substance, the (dis) agreements 
between the derivative expansion and the world line formalism observed in ref. |H] do not 
reflect short-comings of either of these approaches but merely the inappropriateness of the 
fixed flux configuration to check these approaches against each other. 

6.3 Charge Density 

In this subsection, we concentrate on QED 2+1 with a single two-component fermion since for 
a four-component fermion the charge density vanishes identically. 5 Since the charge density 
is a component of a conserved current, it is not renormalized by quantum effects. Therefore, 
we can diagonalize the Dirac Hamiltonian within a Hilbert space that contains a large but 
finite number of states. We then simply compute 



5 For two-component spinors we have tr2 (7 M 7 P 7cr) = ±2ie ppcr where the sign depends on the chosen rep- 
resentation. The four-component spinors may be understood as the sum of these two possible representation 
and hence obey tr4 (7 M 7 P 7o-) = 0. 




(46) 
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Figure 11: The charge density computed from eq. for various locations of the return 
flux and strength of the magnetic held. The return flux is at 6/m in the left panel and at 
8/m in the right panel. 

where ej and ^fj(x) are the eigenvalues and eigenspinors respectively. The technical details 
of this calculation are summarized in Appendix C. We have verified the stability of the 
sum (I46j) with respect to variation of the cut-off that restricts the Hilbert space. In figure ITT1 
we display the resulting densities for the background field defined in eq. (33). Clearly the 
charge density is localized at the regions of the flux tube and return flux. As these regions 
are separated, so are the peaks of the charge density. While the total charge vanishes, the 
integral up to some intermediate point r M « R/2, with p{r M ) « gives 



exactly as expected from a single vortex located at r = [HIE]. In figure we display 
the charge density obtained from the same calculation for background fields without return 
flux. 6 We observe that in the vicinity of the central vortex, the resulting densities with and 
without return flux are identical, again verifying that we may consider the vortex as part of 
a configuration with vanishing total flux. 

It is illuminating to study the charge density using the Green's function formalism as 
we did for the energy density in the previous Section, using the results of Appendix A. For 
background configurations with vanishing total flux, the solutions to Dirac equation come 

6 The calculation as outlined in Appendix A is a unitary transformation on states describing a CP- 
invariant spectrum. Thus the computed total charge will always vanish. However, in the case with non-zero 
net flux, the compensating contribution arises from a peak at a position that in the numerical treatment 
corresponds to spatial infinity, as shown in figure IT51 and is thus considered unphysical. 
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(47) 
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Figure 12: The charge density computed from eq. iBfjj) with and without return flux. In 
both cases eB(0) = 12m 2 . 

in pairs with opposite signs of the energy eigenvalues. From equations (|56|) and (|58p in 
Appendix A we then find that the scattering states obey 




for arbitrarily small but non-vanishing 5. Hence they do not contribute to the charge den- 
sity. Of course, this argument applies to any order of the Born expansion as well. However, 
to apply our Green's function formalism we have to subtract enough Born terms that inte- 
grations in the complex momentum plane become well defined. For the charge density in 
D = 2 + 1 only the first-order subtraction is required. The corresponding Feynman diagram 
gives 

P CO 

PfdO) =rT \ dp 
Jo 

Since we argued that the scattering Green's function does not contribute to the charge 
density we expect the charge to be given exactly by this Feynman diagram. This result 
is confirmed by numerical calculations, shown in figure ED We also compare to the local 
contribution, , which is the leading term of the derivative expansion to PfdM- We see 

that the charge density contains sizable non-localities, mostly in the vicinity of the return 
flux where the background is not slowly varying. 



Jo 



~dr r ^ ' 



arctan ( Jo(p r ) 
\2m/ 



(49) 
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Figure 13: The charge density computed from the Feynman diagram (149)) compared to the 
calculation from eq. (146]) and the local contribution ■ Here eB(0) = 12m 2 and the 

return flux is at 8/m. 

7 Summary 

We have studied exact one-loop quantum energies and energy densities of static electro- 
magnetic flux tubes in three and four spacetime dimensions. To this order only fermion 
fluctuations contribute to the vacuum polarization energy and energy density, which we 
compute exactly from scattering data. In general, this quantum contribution contains ultra- 
violet divergences and an important feature of our approach is that it allows us to impose 
the standard renormalization conditions of perturbative quantum electrodynamics. Even 
though the calculation in three spacetime dimensions does not suffer from such divergences, 
a meaningful comparison between three and four dimensions can only be made when identi- 
cal renormalization conditions are imposed. Thus we must include a finite renormalization 
of the three-dimensional result. Furthermore, the use of scattering data to compute the 
vacuum polarization energy of an individual flux tube leads to subtleties arising from the 
long-range potential induced by the flux tube background, which does not satisfy the stan- 
dard conditions of scattering theory. Consequently, the scattering data do not necessarily 
have the standard analytic properties and the phase shifts are discontinuous at small mo- 
menta. We have circumvented these problems by considering field configurations in which 
the flux tube is embedded with a well-separated return flux so that the total flux vanishes. 
We have constructed a limiting procedure in which this return flux does not contribute to 
the energy, enabling us to compute the energy of an isolated flux tube. While the return 
flux can give a nontrivial contribution to the charge and energy densities (which integrates 
to zero in the case of the energy density), such contributions are well separated from the 
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flux tube and thus easily identified. Thus we have an unambiguous definition of the energy 
and charge of the isolated flux tube even in the embedded configuration. This embedding 
is analogous to considering a kink-antikink rather than an isolated kink to avoid boundary 
effects in one dimensional 4 theory [To] . 

We do not see qualitative differences between three and four dimensions for either the 
energy or energy density, once identical renormalization conditions have been imposed. How- 
ever, we stress that renormalization in the case of three dimensions proved essential to this 
result because the (finite) counterterm contribution turned out to be large, thus causing 
sizable cancellations in the final result. 

We have tested various approximation schemes to the vacuum polarization energy of the 
flux. Though we have observed convergence of the perturbative expansion for small fluxes, 
we have concentrated on the derivative expansion. While we find good agreement between 
our exact calculation and the derivative expansion approximation for the total energy, there 
are discrepancies for the energy density, which arise because the derivative expansion result is 
an approximation to the effective action, which differs from the appropriate matrix element 
of the energy momentum tensor by a total derivative. 

This study gives an initial step toward understanding flux tubes and vortices in more 
complicated theories. The next step would be to consider the complete U(l) model, including 
a Higgs field with spontaneous symmetry breaking, which would describe flux tubes in a 
Type II superconductor. It could then be generalized to non-Abelian gauge bosons for the 
study of Z-strings in the Standard Model, which are flux tubes in the Z gauge boson field. 
The methods developed here could then be used to determine whether quantum corrections 
stabilize the classically unstable Z-string. 
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A Energy Density from Green's Function 

In this Appendix we describe the construction of the exact fermion Green's function in the 
flux tube background in terms of scattering data. We then use it to compute the quantum 
energy density due to the fermion fluctuations. 
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We will begin by considering four-component Dirac spinors \I/ in D = 3 + 1 dimensions. 
The calculation for two-component Dirac spinors in D = 2 + 1 dimensions is then a simple 
manipulation. For the electromagnetic field = (0, A v e^)^, the 4x4 Dirac equation 
decomposes into 2x2 blocks 



m Hi 
H 2 —m 



(50) 



where we now use the Bjorken-Drell representation for the Dirac matrices. Since the system 
is translationally invariant in the z direction, we can choose \l/ to be an eigenstate of z 
momentum with eigenvalue p. Then the 2x2 block becomes 



H 2 
L 



p Lt 
L —p 



-le 



d r + eAJr) + ~cL 



ic 



(51) 



where r = ^ x 2 + y 2 . Accordingly we decompose the four-spinor into two two-spinors \I> 
(£ m , which obey the second-order equations 



y + isl o 





Hfr = (ui-m) H£ = (u, 2 - m 2 ) n with H> = ( 
and similarly for £. Thus we only have to diagonalize L^L and LL^ by solving 
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k 



where h(r) = erA v (r) = jF/(r), as in eqs. (0) and (J3J. Note that £ takes positive and 
negative integer values. 

The equations in ()53|) represent problems in ordinary scattering theory, with potentials 



^(r) = ^-2(/ + l)^-^ 



(54) 
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and their solutions possess the standard analytic properties in the complex k plane if these 
potentials obey the usual conditions [T5]. This is the case for the embedded flux tube, where 
we can use the techniques of ref. j23] to construct solutions for gf' that obey the Jost 
and scattering boundary conditions. Since the rotation to the imaginary k axis removes 
the explicit contribution of any bound states to the energy density [S], it is sufficient to 
concentrate on the scattering solutions. 

We construct the scattering eigenspinors of eq. (joTi|) in analogy to the plane wave solution 
to a free Dirac field. The eigenvalues of the full Dirac problem are uj = ±-^/ 'k 2 + p 2 + m 2 , 
and we define the positive quantity E = \Jk 2 + p 2 + m 2 . Spinors with eigenvalue u = +E 
are 
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E+m g^(k;r)e^J 

( gf\k-r) 


-L^gf\k-r)) 


gf\k;r) 
9?\k 



E+m 
( 



ie-tv T (+)J 2 )l 



\ 



Eh9?\k;r) J 



where we have introduced the linear differential operators 

+ 1 - h(r) 



4 +) =d r + 



and 



Li") = d r ~ 



i - h{r) 



(55) 



(56) 



These equations exhibit the symmetry 

g { Xi(k;r)[h]=gf\k;r)[-h]. 
Similarly, the spinors for the anti-fermions with eigenvalues uj = —E are 



/ 



$f\k,p; r, z) 



E + m 
AtxE 



E+^Pl \k;r) 
' l[~\f\k\ r 
gf\k;r) 



E+m 



\ 







(57) 
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E + m 
AttE 



e i(l+l)<p Q ipz 







(58) 



Assuming the normalization condition for boson wavefunctions as in ref. 



rdr gf*{k'; r)gf\k, r) = -5{k' - k) i = 1, 2 



(59) 



the spinors obey 



d 2 r ^' ] \k\p- r)*f\h\ r) = 6 a rJ e ,A5(k' - k) . 



,(«) 
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(60) 



The solutions in the two-dimensional case are obtained by setting p = 0. Then the spinors 
are also eigenstates of diag(l, —1, —1, 1), reflecting the decomposition of the four-spinor into 
two decoupled two-spinors. 

We would like to compute the renormalized radial energy density, 



e(r) = 2Trr(n\T 00 (x)\n) 1 



(61) 



Here \Q) denotes the vacuum of fluctuating fermions that are polarized by the background 
electromagnetic field, A^x) and 



r 00 (x) = V(x) (-iYdi + eA+m) #1 



(62) 



is the energy field operator. The polarized vacuum is annihilated by the operators af\k,p) 
and bf\k,p) that appear in the decomposition of the field operator 



y(x) 



dp 



dk 

JtTe 



oo I 

^[^ s) (A:, P ;r)e-^af(fc, P ) 



-oo s=l 



+y<f +2) (k,p ] r)e lEt b<f ) \k,p)] . (63) 



Again, we have omitted potential bound state contributions. 

The ultraviolet divergences in the matrix element (J61)) are canceled by counterterms 
that we determine by suitable renormalization conditions. These divergences are at most 
quadratic in the background field, A v . The Taylor expansion in A^ for [e(r)] 2 , the radial 
energy density with the first two Born approximations subtracted, starts at cubic order and 
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is free of ultraviolet divergences. We will therefore first compute [e(r)] 2 and in the next 
Appendix add back the subtracted pieces in form of Feynman diagrams combined with the 
counterterms. We obtain [e(r)] 2 using the standard anti-commutation relations for af\k,p) 
and bf\k,p). Since the spinors ^f\k,p;r) obey the Dirac equation this yields 



e(r)] 2 = -r P d<p [ ^ f° — V WF+p 2 + m 2 fe a)t (fc,p; r)^ a \k,p- r 



-OO = 1 



J 2 

(64) 

We note that Li ^g^\k;r) and gf\k;r) satisfy the same linear second order differential 
equation, as do Lr g\ (k;r) and gf\k]r). To relate one solution to the other, we need to 
consider the boundary conditions. The radial functions that enter the Green's functions (fBTIj) 
obey physical scattering boundary conditions j22j. Thus we identify 



h\. ) gf\k;r) = —ikgf\k\r) and L^gf\k; r) — ikgf\k; r) , (65) 

which simplifies the sum in eq. f!64|) considerably. Our main objective is to express the right 
hand side of eq. ()64|) in terms of the Green's functions for gjp (k; r) and gf\k] r), 

7T J (k + le) 1 - q 2 

by noting that the imaginary part of the Green's function at coincident points may be related 
to the wavef unctions, 

\mGf{r^k)= r 1 -\g\ s \k-r)\ 2 . 
Collecting pieces finally yields the energy density 

where we have also used that \mG ( f\r,r } k) is odd in the momentum k. We have succeeded 
in finding the energy density from fermion fluctuations in terms of bosonic Green's functions 
Gf\r,r' , k). For further details on the decomposition of these bosonic Green's functions 
into Jost and regular solutions we refer to ref. j2Hj. In that paper also the computation 
of the corresponding phase shifts 5^ is described. In the present case, 5f\k) and 8f\k) 
actually refer to the same Dirac spinor. They must be equal, which is a consequence of the 
identification (j63j) . 

It is also straightforward to obtain the expansion in the background field, A v . We 
iterate the differential equations starting from free solutions obeying outgoing wave boundary 



conditions, (k; r) = H^ s _ x (kr) [23] . 
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Once we have removed the ultraviolet divergences by subtracting the leading contribu- 
tions in Atp, the integral in the complex k plane over the semi-circle at infinity vanishes 
and we can compute the k integral by contour integration. Then only branch cuts from 
multi-valued functions (such as \Jk 2 + p 2 + m 2 ) and poles due to bound states contribute. 
In our calculations, any contributions from bound state poles will be canceled by the explicit 
contributions of the bound states themselves, so we need not consider them. Thus if f(p, k) 
is an analytic function in k and goes to zero fast enough at infinity that it does not contribute 
to the integral over the semicircle, the replacement E — > E — f(p, k) does not alter the value 
of the integral (|67|). In particular, we may choose f(p,k) = \Jp 2 + m 2 + k 2 j1\/p 2 + m 2 . 
This modification allows us to carry out the p integration first, 



+00 2 



2m 



-00 s=l 



m 



kG<f\r, r, k) 



(68) 



where Uk = Vk 2 + m 2 and we have dropped pieces analytic in k because they vanish by the 
above argument as well. In the upper half plane, In (a;) has a branch cut along the imaginary 
axis, starting at k — im, where it jumps by m. Therefore the energy density is given as an 
integral along the imaginary fc-axis, 



dt 



-00 2 



2tt 

m i=-oo s=l 



(69) 



where k = it. 

The integral of the Green's function is related to the Jost function F^ s \k) by (22] 



2k 



d 



dr[G (s \r,r,k)} 2 = -— In Fi s \k) 



As), 



for ImA; > 



(70) 



which implies that the twice-subtracted total energy obtained by integrating [e(r)] 2 in the 
form of eq. (jHSj) is given by the phase shift formula 



dr [e(r)]r 



2 JO 27T 



00 2 



=00 s=l 
00 2 



J 2 



2tt 2 



e=-oo s=i 



(71) 



where again terms that do not cause branch cuts in the complex plane have been omit- 
ted. Taking into account that the two phase shifts are identical we may define 5 e+ i(k) = 
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5f\k) 



5f\k) 



to write 



J 2vr 2 m 



which is eq. (jTSjl in the absence of bound states. (On the real axis, the bound state con- 
tribution does not cancel, so to find it requires that we restore the contribution from the 
analytic terms in the momentum integrands.) The fourfold degeneracy for each momentum 
k represents the two possible energies ±Uk and the two decoupled spin channels s = 1,2. 

Finally let us consider the D = 2 + 1 case. In that case we drop the p integral in eq. (jfi7j) . 
yielding 

Integrating along the branch cut we find 



77T 



kGf\r, r,k) 



(72) 



\t2D(r 



/OO J. +°° 2 
1 ^ £=-oo s=l 



(73) 



We recover the expression for the total energy, eq. ([17)1 . similarly to the .0 = 3 + 1 case using 
the spatial integral in eq. (|70|). 



A.l Phase shifts for non-zero flux 

When there is no net flux, the potential, eq. (j3jj) approaches zero as r 
convenient to parameterize the radial fermion wavefunctions as 



?\k;r) = A ) ™H®_ 1 (kr) 



9e 



oo and it is 



(74) 



These ansdtze induce non-linear second order differential equations for the complex radial 
functions p\ (k; r). Imposing the boundary condition lim r 
these differential equations yield the phase shifts [T7] 



/3f (k; r) = 0, the solutions to 



5f (k) = - lim Re [[5[ s) (k; r)) . (75) 

As usual, the Born expansion for the phase shifts is obtained by expanding {k; r) in powers 
of the backgroud field. The solution to the iterated differential equations for (3 ™ {k\ r) then 
yields the n-th Born approximation for the phase shifts 6^ n,s \k) = — lim r ^ ^ e f/?j n,s ^(A;; r ) j 
|E|. 



34 



In the case of non-zero flux, i.e. Hindoo h(r) = T 7^ 0, the situation is more complicated. 
The ansdtze (174)) are changed to 

g { l ) {k-r)=^ ) ^Hfl_ 1+r {kr), (76) 

to accommodate the modified asymptotic behavior. The phase shifts are again defined as in 
eq. ([75)1 . but the Born expansion needs to be modified. It is essentially a series in powers of 
J 7 , but in the product ansatz (fTUj) . both factors depend on T . For the computation of the 
Born series it is therefore more useful to decompose as in eq. (I74j) . 

g\ s \k-r) = e^^Hfl_ l {kr), (77) 

and expand the exponent according to P^-(k; r) = T n p^ (k; r) where n labels the order 
in the Born series. Thus the n-th Born approximation for the phase shift is formally the 
same with and without net flux. However, for T 7^ the background potential approaches 
zero only slowly as r — > 00. Then the numerical treatment requires special care because we 
start integrating from large but finite r^, where it is no longer accurate enough to adopt 
the naive boundary conditions f3jj^(k; r^) = 0. Rather we take j2E| 



$It r oo) — r - ) r ^ 

d 



00 



r=r oc k 

^ ) (*;roo) = 2 j ^, l^ S \k-r) = --L-. (75 



2krl 



to integrate the corresponding differential equations from r = towards r = 0. 



A. 2 Known properties of threshold states 

In this Appendix we review results for threshold states in magnetic backgrounds P3EU]. We 
concentrate on a single two-component spinor in 2 + 1 dimensions, obeying the Dirac equation 

in eq. (|51j). We parameterize the two-component spinor as \1/ = (e lfjp g^p (r) , e^ +1 ^gf\ r i 
At the positive energy threshold, uj = +m, we find the solutions 

g f\r) = /exp(-$(r)) and g^(r)=0, (79) 

where 

$( r ) = f dr'^p- = T f dr'^p- . (80) 
Jo r ' Jo r ' 

At the negative energy threshold, uo = —m, we find 

g f\r) = r""* -1 exp($(r)) and gf\r) = Q. (81) 
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In a system with net flux T ^ 0, f(r') approaches unity as r' — > oo, and thus $(r) — > jFlnr 
as r — > oo. For r^Owe demand /(r') ~ r' 2 , and thus lim,,^ 3>(r) — * 0. Normalizability of 
the wave function thus requires: 

= +m : J^>£+1>0, (82) 
= -m : 0>£> T . (83) 

Hence a positive (negative) flux is required to generate threshold states at positive (negative) 
threshold. The same analysis can be performed for 3+1 dimensions, but there we have two 
two-component spinors of opposite chirality. As a result, for a given flux we find equal 
numbers of states at the two thresholds in that case. 

The threshold states have a very close relation to the famous Landau levels in a constant 
magnetic field B , where the total number of threshold states is infinite. For B > 0, we find 
positive energy threshold states 

5 f(r)=rV^ 2 and gf\r) = Q, (84) 

Analogously for B < we find at negative energy threshold states 

g f ( r ) = r - t - 1 e^ Lr ' 2 and gf\r) = Q. (85) 

Infinitely many threshold states emerge because these wavefunctions are normalizable for 
any value I that allows a regular solution at r = 0. These solutions lead to the standard 
Landau density of states, appearing at the appropriate threshold. 

B Feynman Diagrams 

We begin this Appendix by showing formally that for any static electromagnetic background 
field, the total energy obtained from the spatial integral of the vacuum matrix element of 
the energy density operator equals minus the action per unit time, as one would expect. 
We show that this identification does not extend to the energy density, however. Finally, 
we use these results to compute the first two terms in the Feynman series for the flux tube 
background. 

B.l Feynman series for the total energy 

We would like to develop a Feynman series for the total energy in a static background. It 
will contain divergent terms, so we must carry out our calculation using a gauge-invariant 
regulator, such as dimensional regularization. Since we are considering only static back- 
ground fields, the time coordinate is special. Thus we can use dimensional regularization in 
one time dimension and D — 1 space dimensions, where only the latter may be fractional. 
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The energy density is obtained as the vacuum expectation value of the operator defined 
in eq. (|62|) . In order to generate the Feynman series we adopt the functional language, and 
for simplicity we rescale the vector field to contain the coupling constant e. We have 

(f 00 (x)) = -zTr {(ilp - m)~ l 8{x - x) [-ifdi + A+ m]} 
= —i Tr {{ilp — m) -1 5(x — x) [i7°<9o — ilp + m\ } 
= — i Tr {(ilp — m) -1 5(x — x)27°<9o} 

= -iTili^do^Ip-my 1 5(x-x)} , (86) 

where we have omitted independent terms because they are canceled by subtracting the 
zeroth order term of the Born series, the cosmological constant. The rest of the Feynman 
diagrams are obtained by expanding 

(foo(x)) = -zTr {i 7 °<9 (1 - SAy 1 S5(x - x)} (87) 

in powers of A^. Here S = (i^d^ — my 1 is the free Dirac propagator. The n th term in the 
expansion is 



(T m (x)) 7 



tTr {z 7 ° 'd (SA) n S '6 (x-x)} 



(88) 



Since the background field A^ is static, it is useful to introduce frequency states \u) with 
(u;|j4 m |u/) = A^5{u — to'). Then the expectation value becomes 



(Too(x)) n = -i / ^Tr {i°u[S(u)A] n S(u)5(x - S)} . 



19) 



where S(u) = (7°^ — rf-d — m) 1 and Tr' is the trace over all remaining degrees of freedom, 
spatial and discrete. The total energy at order n is 



d u - l x(T 00 (x)) n 



{siuy^is^AT} 

Z7T 

-/ / ^Tr' {S(uy°uS(u)AS(u)AS(u)A. . .} . (90) 



We use the Ward identity J^S(u) = —S(cJ) , y S(u) to write 



du 



271 

duj 
2^ 



Tr'lco 







duj 

Tr'<j S(u)AS(u)AS(u 



+u 



u)A 



AS(u)AS(u)A. 

. . + Ui 

d 



d_ 

du 



du 



S(u) 



+ 
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(91) 



where we have integrated by parts and used the cyclic properties of the trace. We identify 
the first and the third terms of the above equations to eliminate the derivative terms, 



(92) 



The time interval, T, appears because we have re-established the full functional trace. This 
expansion can then be re-summed 



E = Y,E n = ^Trln[l - SA] = ^Trln[^ 



771 



(93) 



where again independent terms have been omitted. We see that order by order the total 
energy can be obtained from the effective action. 

However, this is not the case for the density. The action density would be 



The terms in this expansion are simpler than those in eq 
Dirac propagators at order n, instead of the n + 1 in eq 



ill {In [1 - SA] 5(x - x)} = ^2(A(x)) n 

(94) 

because they only have n 



(A(x)) 



(A(x)) n = -Tr{[SA} n 5(x-x)} 



Upon inserting 1 = -J^-lu and integrating by parts, eq. (J9"4"j) turns into 



n— 1 



k=0 



duo , 



in—k 



2tt 



8(x — x) \S{uj) 



(95) 



These are the same operators as in eq. ()89j) . but in different orders. The differences are 
characterized by the commutator [8{x — x), S] which has the matrix elements 



(k'\[5(x-x),S}\k) =e t{k - k ' yx 



$ — m ft' — m 



(96) 



These matrix elements have expansions in k — k', starting at linear order. That is, in 
coordinate space they will give total derivative contributions to the energy density (thus 
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leading to the same total energy). When expanding the energy density we deal with three 
operators, S, A and 5(x — x), while the total energy contains only the combination S — A. 
Thus the total only involves two operators, which that can always be brought under the 
trace in the desired order. 

The expansion (|9~%jl will always be of the form 

tr j^x) J d% . . . An-i4(<7i) • • • 4(?n-i)e^ (<?1 - ) n(g 1 , . . . , q n _ x ) J 

with some (n — l)-point function 11(91, . . . , q n -i)- Thus it is semi-local — it vanishes every- 
where where vanishes. 



B.2 Feynman diagrams for energy density 

In Appendix[X]we have computed the energy density with the first two terms of the expansion 
in subtracted. This subtracted quantity was then amenable to contour integration in the 
complex momentum plane. Now we have to add back in the terms we have subtracted, which 
we will do in terms of Feynman diagrams, treating renormalization the standard way. 

We start from the formal expansion in eq. ([89)1 . The first-order contribution is given by 



(TooOr)>i = -1 J dkJhtr {70^1 - m)- x A\k x - k 2 ){ty 2 - m)- 1 ^ 2 "^ J 

= -iJ dldqA^q)e- ic >- x J d£ [l 2 - m 2 + f (1 - £) ? 2 ] ~ 2 

xtr { 7o / [/ + ^ + m ] 7 ^ + (^-l)^ + m ]} . (97) 

where dk = for momenta and dx = dPx for coordinates in D dimensions. Here we 

have already used that the background field is static, so that A^(q) oc S(qo), which allows us 
to set q° = inside the integrand and therefore the shift k\ = I + ^q does not affect the time 
component. We observe that 

1° [I 2 -m 2 + C(l - Ofl " 2 = -~ [I 2 ~m 2 + £(1 - 0? 2 ] 1 (98) 
and integrate by parts in l , 



(foo(rc))i = -- j dldqA,{q)e-^ x J d£tr 



I 2 -m 2 + £(l -Qq 2 
-2i I dqq"AJq)e- iq - x I d£ I dl 



1 2^ — 1 



I 2 - m 2 + £(1 - £)g 2 
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= (99) 

because the ^-integral vanishes. As expected from Furry's theorem, the first-order contri- 
bution vanishes. 

The second-order contribution is 

(f 00 (x)} 2 = -zTr {i lo d SASAS5(x - x)} 

= -2i [ dqdpA fM (p)A I/ (q-p)e- iq x I d£ [ drj (100) 



x dl tr 



jo jo 

^°7o(/ + jj + m)^{]t + + m)Y(} + $ - ^ + m) 
[I 2 -m 2 + A] 3 

where we have introduced the abbreviations 

K = & li + vq?, A = Z(i-Op 2 + v(i--v)q 2 -^vp-q = & 2 + vq 2 -k 2 i (101) 

and again used that the background field is static to set p° = q° = inside the integrand. 
We integrate by parts in 1° to obtain 

(T 00 (x)> 2 = - % -j dqdpA ll (p)A l/ (q-p)e-^W(p,q), (102) 

where 

u^(p iq )=f\ f d v [ - — IL—ixIwu + t-t+mwy+t-i + m) 

Jo Jo J [I 2 - m 2 + AJ [ £ 

+7°(/ + ^ + m) 7 V7"(/ + ^-^ + ^)| • (103) 

Obviously, H ,iu {p, q) contains ultraviolet divergences as D — > 4. However, IP^p, q) should 
become finite by merely adding the conventional counterterm proportional to F 2 U ; no other 
counterterms are available in this theory. To see how this result emerges from the equations, 
we evaluate the trace in eq. (|103|) and keep only terms that do not vanish under the /-integral 



tr{ . . . | = trJ 2 (^j^ 2 + m 2 -A S j + 2A7V 



+27*0 - fly (jf + -fW-f-fW (104) 
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and the quadratic divergences drop out, giving 



n""(p, q) = f di f~" drj I dl ^ iP '^l (105) 
Jo Jo J [I - m + AJ 



with 



^(p,g;e^)=tr{2A7V + 2 7 ^(^-^(^-|) + A 7 V7V+7 ^7V(^-^)). (106) 

To deal with the logarithmic divergence, we first consider IP^Qo, 0). In that case the rj- 
integral becomes trivial and yields a factor 1 — £. The ^-integral simplifies due to the 
identity 



f dze(i - o/(e(i - o) = f dm - e) 2 /(e(i - o) = \ I dm - o/tea - o) > 

Jo Jo z Jo 

which leads to 

dl 



(107) 



u^ip, o) = 8 ( P 2 g^ - pv) deea - o y 



2 • ( 10 ^ 

[P -m 2 + £(1 - Op 2 } 



This has the form of the unrenormalized second order contribution to the total energy, cf. 
eq. ([103)1 . Therefore one might assume that in the decomposition 

(f 00 (x)) 2 = - t - J dqdpA,{p)A v {q-p)e-^{W^p^) + {W v {p ) q)-W v {p^)}} , (109) 

the divergence from the first term in curly brackets would be canceled by the usual coun- 
terterm and the difference II' u '(p, g) — W{p, 0) would give a finite result on its own. The 
latter conjecture can easily be proven incorrect, since the traces in eq. ()1U5|) are different for 
q and q = 0, hence the integrands in the IP^Qo, q) — IP^Qo, 0) do not cancel at large I. 
Nevertheless, let us continue and consider the first term in eq. (|109j) . 

T x {x) = - % -J dqdpA,(p)A u (q-p)e- tc > x U^(p,0) (110) 
= -U J dqdpA^p)A u (q - p)e-^ x (p 2 g^ - p»p») U(p 2 ) , 

where, as usual 

- 1 r dl 
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Note that the piece Ti(x) is semi-local — it vanishes everywhere where A u (x) is zero because 
J dqA u (q — p)e~ iq ' x = e~ ip ' x A u (x) . Furthermore, eq. f|l 1(J|) is the second-order contribution, 
A 2 (x), in eq. CT . 

Next we add the usual local counterterm determined at the renormalization scale M 2 . 
(For the calculations reported in the main text we used on-shell renormalization conditions 
corresponding to M 2 = 0.) This term, 

£ ct = -2ffl(M 2 )F /li/ (x)F"(x) (112) 
= -4T1(M 2 ) {<9 M [A v {x)d^A v {x) - A v {x)d v A»{x)} - A v {x) [g^d 2 - d»d v ] A^x)} , 

does not make the piece T\{x) finite because of the surface term. 7 Rather we find 

T^x) + e ct (x) = -4dA u (x) J dpA^p)e~^ x (p 2 g^ - p»p v ) n R (p 2 , M 2 ) 

+4zn(M 2 ) J dpdqA^p)A v (q - p)^ {q»p v - g^p ■ q) , (113) 

where the renormalized polarization tensor is 

n R ( P 2 , m 2 ) = Ti( P 2 ) - n(M 2 ) . (ii4) 

Since the theory is renormalizable, the divergence in eq. ()113|) must cancel the divergence 
in the difference U tJ,u (p, q) — Tl^ u (p, 0). To evaluate that difference we have to compute the 
numerator trace ([106)1 . Because the background is static and its time component vanishes, 
we may replace ^ Q ^^^°^ u — > —7^7^ and 7°^7 M 7°7 iy (^ — </i) — ► ^7^7" — fi). We obtain 

N^(p, q- e, 77) = 4 [£(3 - 20p 2 + (2 V - 4fr - l)p ■ q + 77(1 - 2/7)g 2 ] 
+16£(£ - 1>V + 8/7(277 - l)qVq u + 8^p u 
+4 [(£ - 1)(2t7 - 1) + 2^ - f] (pV + ?V) • (H5) 



Using 



f dr ] [N^(p,q-,Z,r ] )-N^(p,0-,Z,r ] )} = Uq^-p-qgn , (116) 
Jo Jo 



we have 

•1 r 1 --- 



^ Jo Jo J [t - tttJ 

which indeed cancels the UV divergence in eq. (|113|) . 

7 This result also indicates that the expansion l|94|l cannot be used together with the (standard) countert- 
erm Lagrangian, eq. (|113|l . 
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Putting everything together, we have 



(f 00 (x)) 2 + e ct (x) = 4ffl R (M 2 , 0)«9 M [A v (x)d^A"(x) - A u (x)d v A^x)} 

-ua v (x) J dpA^e-**-* (fg^-v^f) n R (p 2 ,M 2 ) 

"2(^)2 / dqdpA,(p)A u (q-p)e- i ^ * drj 

x { N^(p, q; £, ^)ln^^ - A^(p, 0; ^ln ^^^ 



■ 2ffl R (M 2 , 0)F^F^ - AiA v {x) j dpA^p)^ (p 2 g^ - jSp") Tl R (p 2 , 0) 
■— ^ J dqdpA„(p)A u (q - p)e-^ x j\ jT ^ 

x { g; £, //Jln^-^ - A^(p, 0; £. ^In ™* ~ ^ ~ ^ \. (118) 



We define 

8vrin R (M 2 ,0) = Fi(M) = — - 
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2m 2 


" 2^ 


18 " 


f 3M 2 


-i 
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2m 2 


2^ 


18 ~ 


3p 2 



arcsin 

\2m J 



ln- 



18 3p 2 6p 3 ^4m 2 +p 2 y^ 2 + p 2 — p 



,(119) 



for M < 2m with p = |p | and 



2n 

h(p) = / drJ (pr)h'(r) , 

p 



so that 



A (p) — ih{p)z x p , 
where p is in the x — y plane. We finally find 

2 



(120) 
(121) 



(foo(x)> a + e ct (x) = iF 1 (M)(^) + ^ 



p 3 dpF 2 (p)h(p) Ji (pr) 



■—^5 y qdqJ (qr) J dph{p) ^ 



2tt 



h(\J p 2 + q 2 — 2pq cos0 ) 



\/p 2 + q 2 — 2pq cos0 
x [(p 2 - pqcos(f)) Zi(p, q, cos0) + (pgsin0) 2 Z 2 (p, g, cos0)] . (122) 
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We have used the Lorentz gauge condition, which in momentum space Hows us to substitute 
p u A u (q — p)^ q u A v (q — p) and p^A^{p) — > 0, so that the only tensor structures are 5 mn and 
q m q n - They multiply Feynman parameter integrals 



Zi(p,q,z) 



The rj integrals 



/' 

Jo 



dd f (3 - 20p 2 /o + [(2 - 4£)/i - Io] pqz + {h - 2I 2 )q 2 



Z 2 (p,q,z)= I ^[/o + 4(e-l)/i + 4/ 2 ] 
./o 



i(p,q,z;0 = / 

./0 



d?7 7] n In 



i+e(i-o4+^(i-^)4- 2 ^^ 



(123) 



(124) 



can be performed analytically. We define 7] = q — 2£p,2 and s = ^ Am 2 + 4£(1 — £)p 2 + ??q 
and find 



J = (1-0 ^ 



m 2 + e(i-0(p 2 + r-2p^) 



2 i „2 



770 n rn + £(1 -£)(p +q - 2pqz) s 



-^ln 

2q 



m? + £(1 — £)p 2 



H In 

2g 



m 2 + £(1 — £)G° 2 + <? 2 — 2pg^) 



[2g(l - - % + s] (?7o + s) 



[2q{\ -Q-no-s] (770 - s) 



h = 



Vo + s 2 ^ m 2 + £(1 - £)(p 2 + g 2 - 2pg^) + ^o£ ln 
8q 2 m 2 + £(l-£)p 2 Aq 2 

m 2 + (p 2 + q 2 -2pqz) 2 



[2q(l - ~ Vo + s] (Vo + s) 



[2q(l - - 770 - s] (770 - s) 



^(l-0 3 {ln 

_l (1 _ 25b_i (1 _ i+£! 

6 V W g 6 V W g 2 
77o f^p | ^m 2 + £(1 - 0(P 2 + g 2 - 2pg^) 



+ 



8g 3 
s 

8q 3 



+ 77o In 



m 2 + £(1 — Op 2 
[2g(l-0-77o + s] (770 + s) 



[2q(l-0-Vo-s] (770 -s) 



(125) 



Although these formulas look awkward, they can straightforwardly be included in a nu- 
merical program. The separation of IP^p, 0) was advantageous to integrate over the loop 
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momentum. Once that integral was computed, the ^-integral in eq. (|118|) could be done. It 
then turned out that this term exactly canceled the semi-local contribution. That is, 



(T 00 (x)) 2 + e ct (x) = 2zU R (M 2 , 0)F flu F^ 



2(4vr 
This turns into 



1 r „ r 1 rt-t m 2 - A 

— 2 / dqdpA^MQ - P>- iq X / « / d V N^(p, q; £, V )ln — . (126) 

br) 2 J Jo Jo rn 1 



(Too(x)) 2 + e ct (x) = —F\(M) 

71 



qdqJo(qr) / dph(p) 



2tt 



h(\/p 2 + q 2 — 2pqcos(f 



647T 5 J J Jo a/p 2 + q 2 - 2pq cos0 

[(p 2 ~ VI cos0) Zi(p, q, cos0) + (pq sin0) 2 Z 2 (p, q, cos0)] . (127) 



x 



Here 



Zi(p, q, z) = [ d£ {£(3 - 2i)p 2 h + [(2 - 40/i - h\ pqz + (h ~ 2I 2 )q 2 } 
Jo 



(128) 

while Z 2 is the same as in eq. ()123j) . Though eq. ()127)1 looks simpler than eq. ()122|) it is not 
obvious whether it is easier to handle numerically because in eq. ()122j) the semi-local and the 
total derivative terms have been separated. The latter vanishes under spatial integration, 
and tends to be smaller because in eq. ()123p the leading p 2 lnp terms cancel for Zi, but 
they do not for Z\. The numerical computations presented in the main text are based on 
eq. flT22jl . 

The case D = 2 + 1 is similarly given by eqs. (jl()2j) and ()1()5|) . In particular the numerator 
trace ()115j) is identical, since we chose to work with four-component spinors. However, now 
the loop integral ()105|) is finite and there is no need to split up n^p, q) into n MI/ (p, 0) and 
n MI/ (p, q) — n Mi/ (p, 0). Otherwise, the result is very similar to eq. 



irr ( \\ 1 f a it \ I a u \ ^MvV + <? 2 -2pgcos 
(Too{x)) 2 = — —7 I qdqJo(qr) / dph(p) 1 ' 



327T 4 J J Jo \/p 2 + q 2 — 2pq cos0 

x [(p 2 — pq cos0) Zi(p, q, coscp) + (pq sin(j)) 2 Z 2 (p, q, cos^)] . (129) 

While the integral Z\ is different from the D = 3 + 1 case, splitting up n^Qo, q) does not 
affect Z 2 , 

Zi(p, q,z)= [ d£ {£(3 - 20p 2 /o + [(2 - - h] pqz + (h - 2I 2 )q 2 } 
Jo 
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Z 2 (p,q,z)= [ de[/ + 4(e-l)/i + 4/ 2 ] . 
./o 

Of course, the r\ integrals reflect the reduced dimension 







h = 






h = 


g 2 { 


h = 


-3t7o 


q 3 



ni 



e ct (x) = -F 3 (M)F^ = -2F 3 (M) [ - 



(130) 



I n (p, q, z; £) = / d?7 ^ • (131) 

Jo [m 2 + £(1 — £)p 2 + 77(1 — 7])g 2 — 2^7]pqz] ' 

These integrals can be computed analytically, 

. fr) -2q + 2£q\ . fr) 

arcsm I J — arcsm ^ — 

h = ^ { V™ 2 + e(i - 0(p 2 + q 2 - 2pqz) - V™ 2 + la - Op 2 } + f% , 
2 + 1(1 - 0(p 2 + q 2 - 2pgi) - Vm 2 + |(i - Op 2 } 

+^Vm 2 + £(1 - 0(P 2 + <7 2 - 2pgz) + ^|^- /o ' ( 132 ) 

with 770 and s already defined above eq. (|125j) . 

Finally we allow for a finite counterterm as in eq. (|113j) . This adds 

h'^ 2 



™ 2 - C(l - OM 2 
m 2 ( M ( M 2 \ , 2m- M . 

+ 4 + 7^ 1^-— — - (133) 



32vrM 3 [ m \ m 2 J 2m + M 

to the energy density. This regulator has a finite limit as M — > 0, specifically limA/^o F 3 (M) 
1 

24tt' 



C Charge Density 

In this Appendix we describe our computation of the charge density matrix element 

p = (n\^(x)^(x)\n) (134) 

in the vacuum that is polarized by the flux tube background. Here we consider only the case 
D — 2 + 1. The charge density is finite without renormalization and is not affected by the 
finite counterterm. It can therefore be reliably computed in a simplified approach in which 
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we put the system in a large box and sum the contributions from all modes up to a finite 
momentum cutoff A. Since no renormalization is required for this quantity, saturation of the 
matrix element (j!34|) is observed for large but finite A. We follow the method used to obtain 
the soliton in the regularized Nambu-Jona-Lasinio model [29 (but with two-component 
spinors in this case). We start with the spinors 

*(o) M _ ^ ( u ^(r) \ fl351 



where Ji(qr) are Bessel functions. Discretization is achieved by the boundary condition 
Je(qe n R) = at a large radius R. This boundary condition ensures that no flux runs through 
the boundary at r = R and determines the allowed momenta q^ n . Assume there are Ni such 
momenta below the cutoff. Then we label the energy eigenvalues as 



E in = -\ qL +m 2 , n = l,...,Nt 



Ein = + ™ 2 , n=N t + l,...,2N t , (136) 

where n runs from 1, . . . , 2Ng. The normalization coefficients are easily verified to be 



R\Je+i(qe n R)\\ E in ' ^ ^ 



with the momenta qi n assigned as in eq. fjl36jl . Note that there is no problem with the square 
roots in eqs. (j!35j) and ()137j) because the definition of Eg n ensures that the arguments are 
positive for any value of n. 

These states diagonalize the free Dirac operator in D = 2 + 1 dimensions. The interaction 
associated with the flux tube background is 

h(r) ( e~ i(p 

We now evaluate the matrix elements 



Hint = — ( e i<p o ) • ( 138 ' ) 



(m£'\H + H int \n£) = 5 U > (^E en 5 nm + J drh(r) [u f , m {r-)k n {r) + k m (r)ue n (r)]\ (139) 

which by construction are diagonal in the angular momentum £. Diagonalization of the 
Hamiltonian matrix ()139j) yields eigenvalues and eigenstates = X/n=i ^™^in ■ The 
charge density is then simply given by 

£= — OD fl=l 
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oo 2N e . , 2N e 

= EE E cS^S wo«6.w + momo] . (140) 

£=— oo /x=l n,m=l 

In order to mitigate unphysical boundary effects 8 we subtract the free charge density in 
the same geometry, p^(r) = ^2 en sgn(E in )^l n (r)^£ n (r)/2. The difference p(r) — p(°\r) 
smoothly approaches zero as r — > i?. 

This discretization procedure is not unique. However, a different choice would only cause 
a change in the unphysical boundary effects as R — ► oo. Hence, subtraction of p(°)(r) also 
removes ambiguities in the choice of the discretization condition. 
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